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We study the charge transfer dynamics between a sihcon quantum dot and an indi- 
vidual phosphorous donor using the conduction through the quantum dot as a probe 
for the donor ionization state. We use a sihcon n-MOSFET (metal oxide field effect 
transistor) biased near threshold in the SET regime with two side gates to control 
'^ ■ both the device conductance and the donor charge. Temperature and magnetic field 

^ I independent tunneling time is measured. We measure the statistics of the transfer of 

electrons observed when the ground state D^ of the donor is aligned with the SET 

>■ , states. 
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It has been proposed that the manipulation of electron spins bound to single donors rep- 
resents a solid platform for quantum information processing.^ Building such a system will 
require a profound understanding of the donor electrons properties, such as spin coherence 
and decay times, and hyperfine coupling to the nuclear spin environment. While relevant 
information can be extracted with experimental techniques, such as conventional electrical 
measurements and electron spin resonance (ESR) spectroscopy, sensitive to the donor in the 
bulk of the semiconductor, many properties of interest differ when donors are included into 
nanostructured devices because of the presence of material interfaces,^ confinement in po- 
tential wells, and high electric fields.'^ The effect of these parameters can only be evidenced 
in devices where experiments are performed at the single donor level. Recently there has 
been significant progress in the detection of charge and spin at the single donor level in 
silicon^ on devices achieving capacitive and tunnel coupling between a few locally implanted 
donors and a silicon SET working as a charge detector. We have also studied in the past 
the interaction of single donors in the channel of nanometer scaled silicon MOSFETs (metal 
oxide field effect transistors) with quantum dot states at low temperature.^ Here, we achieve 
the transfer of an electron from the quantum dot to a donor in the vicinity of the channel 
by controlling multiple gates and detect the modification of the charge state by observing 
the current flowing in the transistor. In this work we improve the capability to detect the 
charge state of the donor and achieve real time detection of the charge state by employing 
an ad-hoc developed device design. While based on a working principle similar to what 
proposed in Ref. 6, we focus on the development of a more simple and robust technology 
employing a single level of metallization and a substrate uniformly doped with phosphorous. 
The separate control of multiple gates allow biasing the device in a region where the ioniza- 
tion of only one donor is energetically allowed. From time averaged traces of the capture and 
emission process we measure the average tunneling time that is observed to be independent 
on both temperature and magnetic field. From real time traces of the charge transfer we 
measure the average tunneling rates and the process statistics in the framework of the full 
counting statistics (FCS) following the model introduced in Ref. 7 and developed in Ref. 8. 
The tunneling rates to and from the quantum dot reflect the internal energy level spectrum 
of the quantum dot and its occupation. As predicted, the capture process is measured to be 
sub-poissonian as required by the Coulomb blockade regime that allows the capture of only 
one electron in the donor. 




v^f (V) '■' 

gate ^ ' 



FIG. 1. (color online) a) SEM image of a device nominally identical to the one characterized and 
schematic representation of the circuit connections, b) Source-drain current vs. main gate bias, 
with the side gates grounded, measured near threshold at 4.2 K. 



The experiments were performed on a silicon n-MOSFET with an engineered structure of 
gates. The device was realized with a standard process however, in contrast to what used 
for the realization of silicon quantum dots in other experiments, ^'^'^'^ the functionalities were 
obtained with only one level of metalization. A scanning electron microscope (SEM) image 
of a device nominally equal to the one here reported is shown in Figure 1(a). The main 
gate is shaped to form the W = 40 nm and L = 300 nm channel of the transistor when 
a sufficiently positive gate voltage Vg is applied. Near the thin channel two side gates are 
realized at 60 nm from the main gate. The two side gates are biased at the same potential 
Vside during the experiments. The gate oxide thickness, 35 nm, is chosen to guarantee that 
both main and side gates achieve a good control on the charge in the channel. This thin 
channel transistor is used as a highly sensitive electrometer to detect the charge state of 
nearby donors. The phosphorus donors interacting with the channel of the transistor are 
provided as a uniform doping distribution in the substrate of 1 ^ 5 x 10^^ cm~^ and a total 
of 1 -^ 10 donors are expected to be near the active region of the device. At the low oper- 
ating temperature of the device all donors in the substrate are in the neutral charge state 
(D°) and do not contribute to the source/drain conduction. Preliminary characterization 
was performed in liquid helium at 4.2 K while time transient experiments were performed 
in a ^He refrigerator at 300 mK. The sample was connected to a cryogenic amplifier ther- 
mally anchored to the 1 K-pot of the cryostat. The amplifier is a custom circuit realized in 
standard 3.3 V complementary metal oxide semiconductor (CMOS) technology. With this 
configuration the SET current could be measured with a 34 kHz bandwidth and 2 pA noise 
at full bandwidth. 
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FIG. 2. (color online) Derivative of the source-drain current vs. side gate bias {dId/dVside) as a 
function of main gate and side gates polarization at 4.2 K. The arrows indicate the discontinuity 
lines attributed to donor ionizations. The bold line indicates the region of the donor characterized 
in this work. 



By exploiting the disorder at the interface, causing the formation of isolated potential min- 
ima, when biased near threshold the transistor behaves as a single electron quantum dot as 
evidenced by the drain current (J^) vs. the main gate voltage (V^) trace at 4.2 K reported in 
Figure 1(b). When the electron density in the channel is increased the potential fluctuations 
are screened and the ohmic behavior is recovered. Measuring the Id vs. Vg and Vside (side gate 
bias) and plotting the derivative of the current to the side gate bias, as reported in Figure 2, 
it is shown how both gates control the transistor current. Plotting the derivative of the cur- 
rent evidences the Coulomb blockade peaks that show a fairly regular pattern characterized 
by parallel conductance lines. Along several lines in the d{Id)/d{Vside) vs. {Vg, Vside) diagram 
of Figure 2, evidenced by the arrows, it is possible to observe discontinuities in the otherwise 
regular pattern of the Coulomb blockade peaks. These discontinuities are associated with a 
sudden modification of the potential of the quantum dot caused by the ionization/charging 
of a defect state in the vicinity of the channel, as also reported in Ref. 4 and Ref. 11. The 
lines evidenced in Figure 2 are characterized by different slopes in the Vg, Vside plane and are 
thus associated to distinct defects with different coupling to the two gates. After each cool 
down cycle we observed a shift of the Coulomb blockade features in the transport through 
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FIG. 3. (color online) Detail of the Id vs. Vg and Vgide in the quantum dot region near the ionization 
of one donor acquired at 4.2 K. The solid and dashed lines represent the charge stability points for 
the quantum dot and the donor respectively. Current flow is only observed along the quantum dot 
charge stability line. 



the quantum dot that are attributed to defects at the Si/Si02 interface; however the slope 
and relative position of the features attributed to the donors ionization were reproducible. 
This allowed an unambiguous identification of the donors and, after each cool down cycle, 
is was possible to tune the device to characterize always the same donor. This required an 
adjustment of the bias point. The bold arrow in Figure 2 evidences the feature associated to 
the donor that has been characterized in detail. The features attributed to the charging of 
a single donor could not be observed on similar devices realized on p-tjpe silicon substrates. 
Both types of devices are based on a n-FET and are operated in the same bias range; con- 
sequently in both cases the channel interacts with defects with energy states close to the 
conduction band minimum and would show similar features if those were consequence of the 
trapping at defect states. 

When the device is polarized in the quantum dot regime the change of the ionization state of 
a single donor can cause the complete blockade of the source/drain current. This situation 
is reported in the detail of the 2D gate scan acquired at 4.2 K in Figure 3 where the stability 
lines for the quantum dot and the donor are evidenced. The capacitive coupling between 
the donor and the quantum dot causes the shift of the quantum dot conduction line bringing 



TABLE I. Lever arm factors to the quantum dot and the donor from the main and side gates (in 

eV/V). 





Main Gate 


Side Gate 


Quantum Dot 


0.16 


0.075 


Donor 


0.028 


0.059 



the quantum dot out of the resonance hne. 

The lever arm factors to the quantum dot, (yQD,gate{side) = dfiQD/dVg(side), have been ex- 
tracted near the bias point of Figure 4(a) at the temperature of 300 mK from trans- 
port experiments (not shown here), while the arm factors to the donor, adonor,gate{side) = 
d^donoT / dVg(^side)i Were extracted from the thermal broadening of the ionization transition, 
assuming that the quantum dot has a quasi continuous distribution of states with an occu- 
pation given from the Fermi distribution, as described in Ref. 4. The extracted values are 
reported in Table I and are compatible with the interpretation of the presence of a donor 
that is spatially separated from the quantum dot and more closely coupled to the side gates 
than to the main gate. In all the following experiments, we move the potential along one 
line on the Vg^Vside diagram, specified by the arrow in Figure 4(a), and parallel to one of 
the Coulomb blockade peaks. The origin of the new axis is set so that conduction through 
the quantum dot is only allowed when the donor is ionized: in this condition the average 
current through the quantum dot is proportional to the average occupancy of the donor. 
The new V^; bias is proportional to the energy detuning, AE', between the quantum dot and 
the donor, while moving along this axis keeps the quantum dot chemical potential fixed. 
The lever arm factor to the donor is aQD,x = AE/V^ = (0.009 ± 0.001) eV/V. 
The average tunneling times to and from the donor at large bias detuning (|Aii^| ^ ksT) 
are measured by switching the bias to allow the capture and release of one electron and 
measuring the current through the quantum dot. Averaging several acquisitions the expo- 
nential decay of the current, and thus of the donor occupancy, is reconstructed, as reported 
in Figure 4(b). The decay time extracted from the single exponential fit of the current is 
reported in Figure 4(c) and Figure 4(d), and it results to be largely independent on both 
temperature and magnetic field, proving that transfer of the electron is through elastic tun- 
neling between states with the same spin projection along the magnetic field direction. 
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FIG. 4. (color online) a) Detail of the scan of the source-drain current as a function of the two gates 
polarization at 300 niK. The solid arrow represents the line along which the potential is moved for 
the ID scans (the arrow indicates the direction of V^^ positive, b) Average of the transient quantum 
dot current acquired after the switching of the bias to allow tunneling of one electron into the donor 
state. c,d) Dependence of the average tunneling time as a function of the temperature (c) and the 
magnetic field (d). 



When the detuning energy AE is comparable with ksT electrons are continuously exchanged 
between the quantum dot and the donor. The quantum dot current shows random telegraph 
signal (RTS)^^ reflecting the instantaneous occupation of the donor.. Typical traces and the 
reconstructed signal are reported in Figure 5(a). The average time to enter (exit) the donor is 
related to the tunneling time and the density of available occupied(empty) states for tunnel- 
ing from(in) the quantum dot. Assuming that the quantum dot presents a quasi continuous 
density of states occupied according to a Fermi distribution, the tunneling times can be 

written as in Ref. 13: 

1 



r 



out/in 



(1) 



1 + exp{±AE/kBT) ' 

By sweeping the V^ control voltage, we vary AE and the relative occupancies of the empty 
(high current) and occupied (current blocked) states. The experimental tunneling rates, 
together with a fit according to Eq. 1 are reported in Figure 5(b). The value for Fq cor- 
responds to the tunneling rate that has been measured with the time averaged technique. 
Fq = 1.5 X 10^ s~^. The fit of the data in Figure 5(b) reveals a temperature of the elec- 
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FIG. 5. (color online) a) Typical RTS trace acquired near the ionization point of the donor and 
signal reconstructed via a threshold algorithm, b) Average values for the tunneling rates Tin and 
Tout as a function of the detuning potential V^ and fitting according to Eq. 1. The arrows evidence 
the fluctuation in the Tout that can be attributed to resonant tunneling to the orbitals in the 
quantum dot. c,d) Average number of tunneling events per unit time (c) and Fano factor (d) 
(squares) extracted from the RTS traces. The solid lines represent the theoretical value extracted 
from the FCS model using the parameters obtained from the fitting of the data in (b). 

tron system of T = 1.2 K, higher than the refrigerator temperature of 300 mK. Because of 
the weak electron-phonon coupling the electron temperature is generally different from the 
nominal temperature of the ^He bath and of the crystal: this has been observed both in our 
system on different devices^^'^^ as well as in other systems.^ The oscillation of the value of 
the tunnel rates, evidenced by the arrows in Figure 5(b), can be attributed to the discrete 
spectrum of energy levels in the quantum dot: from the period of the oscillation and the 
previously determined lever arm factor the energy separation can be roughly estimated to be 
300 fieV. In a harmonic potential well AE = h ?7ij~^/~^ where ml is the effective transverse 
electron mass and / is the wavefunction extension. From this equation we obtain an effective 
dot size of 36 nm. Because of the confining potential of the dot is induced by charge defect 
at the Si/Si02 interface this number should be compared with the expected defect density. 
A 36 nm linear distance between defects would correspond to an areal density of 8 x 10^° 
cm~^ which is compatible with the technology used. 
From the real time traces we also extract the first and second order statistics of the tunnel- 
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ing process. The analysis, in the framework of the full counting statistics, is performed by 
measuring the probability function of finding n transitions in a given time to and calculating 
the first and second order moments as in Ref. 13: 

rui =< n >t (2) 

m2 =< {n- < n >Y >t (3) 

where <...>t represents the average over time. The first moment mi of the probability 
distribution corresponds to the average number of transitions per unit time, while the second 
moment 7712 is the variance of the probability distribution. In a perfectly Poissonian process, 
with no correlation between the occurrence of two subsequent events, the Fano factor F = 
m2/mi, is always unitary. This value can be modified by the occurrence of some degree of 
correlation. One notable case is that of electrons tunneling in and out of quantum dots in 
which, due to the Coulomb blockade, only one electron at a time can be confined.^ In such 
systems the Fano factor depends on AE as:^^ 

^=(l + /?^)/2 (4) 

where /3 = 2/ {AE/kT) — 1 and / (x) is the Fermi-Dirac distribution. This value reaches its 
minimum at AE = 0, where Te = Tc and F = m2/mi = 0.5. In Figure 5c and 5d we report 
the values for nii and the Fano factor extracted from the reconstructed RTS traces together 
with the values calculated from Eq. 4 with the same parameters extracted from the fit of 
the tunneling rates. 

Concluding, we have analyzed the dynamic charge transfer between a quantum dot and a 
single donor atom in a silicon device. The real time occupancy of the donor was extracted 
from the current flowing through the quantum dot. The device, consisting in a silicon n- 
MOSFET, allowed the real time measurement of the charge state of neighboring donors. 
We studied the interaction with one of the donors with tunneling time in the order of 
1 X 10"^ s. The tunneling rate is found to be independent on the temperature, indicating 
elastic tunneling processes. RTS traces were measured near the energy equilibrium point, 
where electrons are continuously exchanged between the quantum dot and the donor. The 
tunneling rates to and from the quantum dot reflect the thermally broadened occupation of 
the energy levels in the quantum dot. As predicted the capture process is measured to be 
sub-poissonian as required by the Coulomb blockade regime that allows the capture of only 
one electron in the donor. 
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